Theory of Directed Polymers 



Randall D. Kamien, Pierre Le Doussal^ and David R. Nelson 

Lyman Laboratory of Physics 
Harvard University 
Cambridge, MA 02138 

We develop a theory of polymers in a nematic solvent by exploiting an analogy with 
two-dimensional quantum bosons at zero temperature. We argue that the theory should 
also describe polymers in an isotropic solvent. The dense phase is analyzed in a Bogoliubov- 
like approximation, which assumes a broken symmetry in the phase of the boson order pa- 
rameter. We find a stiffening of the longitudinal fluctuations of the nematic field, calculate 
the density-density correlation function, and extend the analysis to the case of ferro- and 
electrorheological fluids. The boson formalism is used to derive a simple hydrodynamic 
theory which is indistinguishable from the corresponding theory of polymer nematics in an 
isotropic solvent at long wavelengths. We also use hydrodynamics to discuss the physical 
meaning of the boson order parameter. A renormalization group treatment in the dilute 
limit shows that logarithmic corrections to polymer wandering, predicted by de Gennes, are 
unaffected by interpolymer interactions. A continuously variable Flory exponent appears 
for polymers embedded in a two-dimensional nematic solvent. We include free polymer 
ends and hairpin configurations in the theory and show that hairpins are described by an 
Ising-like symmetry-breaking term in the boson field theory. 
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1. Introduction 

1.1. Overview 

The statistical mechanics of directed, interacting lines has received renewed attention 
recently. This problem is, for example, directly relevant to the behavior of high-T c super- 
conductors in a magnetic field |]||. In these systems, above a critical external field H c i, 
the magnetic field penetrates the material in the form of lines, each carrying one quantum 
of magnetic flux. The flux lines can be viewed as "polymers" aligned with the direction 
of the external field up to thermal fluctuations. Due to their mutual repulsion, these lines 
can form various states such as a triangular solid @, or isotropic and hexatic entangled 
fluids Hill • Glassy states are also possible, induced either by local disorder |5j|| or simply 
by very long disentanglement relaxation times |l|J^]. Detailed calculations for flux liquids 
are possible by exploiting a mapping onto the statistical mechanics of boson worldlines in 
two spatial and one timelike dimension [1]]. 

Many other physical systems consist of extended one-dimensional objects aligned in 
one direction [j^PJICfl . Stiff biological macromolecules such as DNA [JTJ , helical synthetic 



polypeptides such as Poly (7 -benzyl glutamate) (PBG) [p!2j p!3|, discotic liquid crystals com 



posed of stacks of disk-shaped molecules [^,|T5[, and micelles of amphiphilic molecules JT6 | 
can all form crystalline columnar phases with in-plane order, as well as nematics with fluid- 
like in-plane order. Although stiffer chains align more easily, some nematic polymer liquid 
crystals can also be formed with chains of relatively low rigidity, by alternating a nemato- 
genic unit with a flexible hydrocarbon spacer fll7] - |l9|| . The transition from isotropic melt 



to nematic is achieved experimentally by lowering the temperature J2JJ or more frequently 
by increasing the concentration. Steric repulsion is sufficient to produce alignment at high 
enough concentration, although many other interactions can be present depending on the 
material: Van der Waals attraction, electrostatic forces pl|,p2|, hydration forces, etc.. Fer- 



rofluids [[23 1 and electrorheo logical fluids [24| are also composed of chains of particles, in 



this case aligned by external magnetic or electric fields. 



In a recent paper [25], X. Ao, X. Wen and R.B. Meyer have presented X-ray scattering 
data on PBG which is in many ways strikingly similar to predictions for neutron diffrac- 
tion by flux lines. These authors stress an analogy between polymer configurations and a 
fictitious "dynamics" of two-dimensional particles moving along the direction of alignment, 
and show that the distinctive "bowtie" scattering contours change character in the limit 
of small momentum transfers. 
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There are, however, significant differences between these directed polymer-like objects 
and flux lines. Although oriented nematic polymers in a solvent wander along a preferred 
axis, just as thermally excited flux lines do, the average polymer direction represents 
a spontaneous, rather than externally imposed, broken symmetry. Even if the monomer 
chains are aligned by an external electric or magnetic field, the lines are typically of variable 
length and need not span the system, unlike flux lines. Polymer nematics can, moreover, 
make relatively low energy hairpin turns, because of the symmetry of the director field 
under n — > —n. Such "backtracking" can usually be ignored for flux lines 

In a recent paper [^6| two of us adapted methods developed for flux lines to these 
systems, taking the above differences into account. We showed that the boson field theory 
of becomes applicable to directed polymer melts upon adding a source term, in analogy 
with the des Cloiseaux trick [jT[ for isotropic polymer solutions. The spontaneously broken 
symmetry of polymer nematics requires, in addition, coupling the boson order parameter 
to a massless fluctuating background director field. 

In this work, we describe these calculations in detail, and discuss as well their va- 
lidity in the dilute limit. Calculations for dilute directed polymers are a straightforward 
extension of results for flux lines JJ] in the case of electrorheological and ferro-fluids. The 
dilute limit is more interesting, however, when the average polymer direction represents a 
spontaneously broken symmetry. Polymer nematics in an isotropic solvent will, of course, 
eventually crumple into an isotropic phase upon dilution J27J. An analysis of the tran- 
sition to an isotropic polymer melt would take us beyond the scope of this paper. We 
shall, however, consider a dilute collection of polymers aligned by a solvent which is itself 
a short chain nematic. As pointed out by de Gennes |§, the Goldstone modes associated 
with the the nematic solvent already lead to logarithmic anomalies in the wandering of one 
isolated polymer. Additional logarithms appear when interactions are included, and a full 
renormalization group analysis is required to sort out the details. Yet another complication 
appears when we allow for hairpins in the polymers as they meander through the nematic 
solvent. We find in this case that an Ising-like symmetry- breaking term appears in the 
boson field theory. 

Table I contains a summary ofthe different types of directed polymeric systems con- 
sidered in this paper. Note that flux lines in high temperature superconductors are the 
simplest case, because (1) they contain no free ends and (2) hairpins are highly disfavored 
by the external magnetic field. Although hairpins are disfavored by the magnetic and 
electric fields necessary to produce aligned chains in ferro- and electrorheological fluids, 
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free ends are, of course, unavoidable. The behavior of polymer nematics and polymers in a 
nematic solvent is complicated both by the existence of free ends and hairpins, and because 
the alignment can be produced by a "soft" broken symmetry instead of an external field. 

We consider here only Gaussian fluctuations about the state which describes directed 
polymer melts in the dense hydrodynamic limit. While this paper was in preparation, 
we learned of interesting work by Toner |28| who introduces nonlinearities directly into 
the hydrodynamic theory of polymer nematics. Toner concludes that these nonlinear 
terms eventually trigger a breakdown of hydrodynamics (i.e. a singular dependence of 
hydrodynamic parameters on wavevector) at sufficiently long wavelengths. It would be 
interesting to see if such a breakdown also occurred in the more microscopic boson theory 
developed here. 

1.2. Model 

We shall concentrate on polymer nematics, regarded as directed polymers interacting 
with a background nematic field. By imposing a magnetic field H, or taking the limit of 
very large Frank constants, we can, if desired, recover results for directed polymer melts 
with an externally imposed direction. We start with a nematic free energy 

Fn = lJ d 2 r ± J dz [^(V^Sn) 2 + K 2 (V ± x5n) 2 + K 3 (d z Sn) 2 + H(5n) 2 } , (1.1) 

where the {K{\ are the usual Frank constants for splay, twist and bend, and 5n(r) = 
[5n x (r ± , z), 5n y (r ± , z)} is a vector representing a small deviation of director field n(r) from 
its average orientation along the z-axis, n(r) w (<5n, 1). We neglect for now polymer free 
ends and hairpin turns, and describe the position of the ith polymer as it traverses the 
nematic medium by a function Rj(z) = (r 7 - (z) . z). The iV polymer lines interact with each 
other and the nematic field via a free energy adapted from || 

(1-2) 

+ litjdzV(\r i (z)-r j (z)\) 

Here k is the polymer bending rigidity, while g controls the coupling between the local poly- 
mer direction and the nematic matrix. This coupling is the only one allowed by rotational 
invariance, to lowest order in dvjjdz and 5n. The potential V(r) represents short range, 
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excluded volume effects, and can be approximated by V(r) ~ Vo5 2 (r). The probability of 
a particular field configuration is proportional to exp(— F/k B T), with F = F n + F p , and 
averages are calculated by integrating over both 5n(r) and polymer configurations {vj(z)}. 

The simplest physical interpretation of the free energy F is of polymers aligned by a 
nematic solvent with Frank constants {Ki}. We believe, however, that F also describes 
dense nematic polymers in an isotropic solvent. As illustrated in Figure 1, n(r) then 
represents a coarse-grained nematic field obtained by averaging over the polymer tangents 
in a hydrodynamic averaging volume. Deviations of the orientation of any individual 
polymer from this average direction are described by the coupling g. This interpretation of 
F is especially appropriate for polymers made of nematic molecules connected by flexible 



hydrocarbon spacers |q,|19| as in Figure 2. In this case 5n(f) describes fluctuations in 
the orientations of individual nematogens, while Tj(z) describes how the nematogens are 
threaded together by the hydrocarbon spacers. The bare Frank constants in ( [1 . 1|) are 
then approximately those of the nematic phase of the unpolymerized nematogens. In 
this picture, we have | r = K^/po, where po is the areal density of polymers cutting a 
constant-z cross section. Note that the potential V(r) represents a scalar interpolymer 
interaction within a constant- z plane. The coupling g contains both a scalar interaction, 
due to the longitudinal modes of the nematic, and a vectorial interaction, due to the 
transverse nematic modes (see Section 6). In going from the bare polymer system to this 
coarse-grained system, we could first match the strength of the vectorial interactions with 
g, and then adjust V so as to get the correct strength of the scalar interaction. 

Our assumption that polymers interacting with "nematic background" field are equiv- 
alent to dense polymer nematics in an isotropic solvent is supported by the hydrodynamic 
approach to correlation functions discussed in Section 5, which gives identical results for 
these two systems in the limit of long wavelengths. 

The second basic assumption underlying our model calculations is that we can neglect 
the term proportional to k in ( |1.2j ). To justify this, note first that the initial two terms 



of ( |1.2| ) define a length || A e a/ ft/g, which for isotropic solvents we identify with the 
"deflection length" discussed by Odijk |H| . The deflection length is the distance a polymer 
wanders along z before it feels the confining effect of its neighbors. In order that the 
polymers order nematically, this length must be less than the polymer persistence length 
£ p = n/k B T. On scales larger than A, the coupling to the background nematic field 
dominates the bending rigidity and we are justified in neglecting the first term in ( |1 . 2|) . 



If the deflection length is known, we can express the coupling g in terms of experimental 
parameters as 



9= TT2- ( L3 ) 



PoX 2 

The equivalence between polymer nematics and polymers in a nematic solvent does 
not extend to the crumpling transition to a more isotropic phase which should occur in the 
dilute limit for polymer nematics. The nematic solvent acts like an ordering magnetic field 
in the latter case, so that the polymer exhibits nematic order at arbitrary dilution. We 
should also mention that refinements in the model defined by ( |1 . 1| ) and (|1.2| ) are required 
to handle polymer free ends and hairpins. These are most easily treated after rewriting 
the theory in a second quantized "boson" formalism, as discussed in Section 3. 

1.3. Outline 

In Section 2 we review the wandering of a single polymer in a nematic solvent, as- 
suming at first that hairpin turns can be neglected. In an external field, the polymer just 
executes a gaussian random walk in the xy-plane as it meanders down the z-axis. When 
the field is turned off, the Goldstone modes of the nematic matrix induce logarithmically 
divergent "superdiffusive" behavior, as first pointed out by de Gennes ||. We then gen- 
eralize the model to include hairpins, using an effective field theory method introduced 
by Cardy The hairpins cause the polymer to crumple and execute an anisotropic 

three dimensional Gaussian random walk. We also discuss polymers wandering in a two 
dimensional nematic solvent. Below the Kosterlitz-Thouless transition of the solvent, we 
find that the polymer exhibits a continuously variable Flory exponent, which is simply 
related to the decay of order in the orientations of the solvent molecules. 

In Section 3, we show how to solve the model when many interacting aligned polymers 
are present via a mapping onto the quantum mechanics of two-dimensional bosons. The 
constraints imposed on the theory by rotational invariance are discussed in Appendix A. 
Adding a source to the boson field theory allows us to obtain analogous results to polymers 
of finite length. In Section 4, we calculate the polymer density correlation functions and 
discuss the renormalized wave vector dependent elastic constants. 



In Section 5, we generalize the hydro dynamic approach of de Gennes [30] and of 



Selinger and Bruinsma [^TJ to allow for polymer heads and tails and show that the results 
agree with the long wavelength limit of our more microscopic calculations. In Appendix 
B, we show explicitly how to derive the hydro dynamic theory directly from the boson 
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formalism. We also use the hydrodynamic theory to discuss the physical meaning of the 
"boson" order parameter used in Sections 3, and to calculate the elastic energy of a chain 
end. 

The behavior in the dilute limit is discussed in Section 6. We construct renormalization 
group recursion relations and show how polymer wandering is affected by both nematic 
Goldstone modes and interpolymer interactions. In section 7 we then introduce hairpins 
and show that an Ising-like phase transition then describes the dilute limit. 



2. A single chain in a nematic matrix 

2.1. Three Dimensions: the de Gennes Approximation 

With the problem of the nematic polymer in mind, de Gennes introduced the following 
free energy for a single chain without hairpins in a nematic solvent in the one Frank constant 
approximation (K = K\ = K2 = K3) [§, 



1 f , fd 2 r 

Fi = r dz [d^ 



+ ^ / (h 



dr 

dz 



<5n(r(z), z) 



+ \k I d s r |Vcm| 2 



(2.1) 



The form of the polymer-nematic coupling comes from the small tipping angle expansion 

21 



of \g 



l-(^.fi(7?U; 



rotational invariance and the discrete n — > —n symmetry. We use 



which is the only leading order coupling consistent with 

~ ( dr 1 W f] 1 \dr\2 
d.z ~ v dz ' ^ll \ L "T" I dz I ' 



and n = (<5n, 1)/ y/l + \5n\ 2 to find that 



dR 

dz 



n 



dr 

dz 



5n 



(2.2) 



which leads to the coupling displayed in (|2.1|). 

For a fixed configuration of the polymer, the nematic matrix is distorted and some 
elastic energy results (see Figure 3). We wish to compute the effective free energy of the 
chain resulting from integrating out the nematic field in (|2.1| ). Although all terms in ([2.1] ) 
are quadratic, this problem is, in fact, quite nonlinear due to the appearance of r{z) in 
the argument of Sn(r(z) : z) in the coupling term. A natural approximation, implicit in de 
Gennes' discussion [|], is to set r(z) ~ 0, and to discuss only the effect of the fluctuations 
of the tipping angle 4j. This amounts to replacing the ^-coupling in (|2.1|) , by: 



dz 



dr 

dz 



dz 



dr 

dz 



-8n(0,z) 



(2.3) 
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where 5n(0, z) is the value of the distortion of the nematic field on the line r = 0. One 
expects this approximation to be better for large K. A perturbation expansion in K~ x of 



the original problem (|2.1| ) shows that for d < 3, divergent integrals appear, and one must 
use the renormalization group to obtain correct results. We defer a systematic treatment 
of this problem to Section 6 and discuss here only the de Gennes approximation, which 
can be expected to be qualitatively correct in d = 3. 

The fastest way to compute the effective free energy of the polymer is to note that 
modified by ( |2.3| ) can be simplified by the change of variable (with unit Jacobian) 
r](z) = dr/dz — 5n(0 : z). Upon neglecting the bending energy term, we find that r(z) is 
the solution of 

^ = 5n(0,z) + r](z) , (2.4) 

where rj(z) is a Gaussian white noise of variance k B T/g. The wandering of the chain is 
then given by integrating ( |2.4| ) over z, squaring and averaging over 5n(0, z), 

(|r(L)-r(0)| 2 ) = ^L+ / / dzdz' (5n(0, z) ■ 8n(0,z')) Fn . (2.5) 
9 Jo Jo 

Due to the Goldstone mode associated to the rotational invariance of the nematic free 
energy F n , the director correlation function decreases like l/\z — z'\ at large separations. 
Upon substituting the full Frank energy ( |1 . 1| ) for the one-Frank-constant term in ( |2.1|) we 
find (upon setting H = 0), 



(tfn(0,*).tfn(0,0)) = *2LE 



1 1 



The long range correlations in the medium then implies via (|2.5|) "hyperdiffusion" of the 
chain in the transverse direction: 

/i (t\ fn\\2\ 2k B T k T3 T(K 1 + K 2 ) 

(\r(L) -r(0) ) ~ L+ L In (L/a) , (2.7) 

where a is a microscopic cutoff. The effective "diffusion constant" (defined by (|r(L) — 
r(0)| 2 ) = 4D(L)L) is finite for finite L, but diverges logarithmically as L — ► oo. 

For K\=K.2=K , this is the result obtained by de Gennes by explicitly carrying out the 
integration over the nematic field H. The renormalization group treatment of Section 6 



(which includes interpolymer interactions) leads to the same result with, however, a factor 
of two difference in the coefficient of the logarithm. 

The nematic solvent is not quite able to produce a finite renormalized diffusion con- 
stant for infinitely long polymers. The polymer is still aligned with the z-axis on large 



length scales, however, since y(|r(L) — r(0)| 2 ) <C L. It is easy to show that the bend- 
ing energy, which has been neglected in these calculations, is irrelevant for a single chain 
without hairpins. 

It is interesting to note that, despite its annealed character, this problem is very 
similar to the problem of random walks in quenched random disorder for which long range 
correlations are known to modify diffusion ||32| . 

2.2. Two Dimensions 

It is interesting to ask what the configurations r(s) are of a chain in a two dimensional 
nematic. Although more difficult than its three-dimensional counterpart, an experiment 
on a two-dimensional surface might be possible: Imagine a long polymer chain with N 
monomers adsorbed at an air water interface which is also covered with a tilted monolayer 
Langmuir-Blodgett film. The projection of the tilted hydrocarbon chains on the the plane 
of the interface plays the role of a director ft, without the inversion symmetry n — > —n. 
We parameterize the director by ft = (cos <f>, sin <p) . Thermal fluctuations destroy the tilt 



order above T c via a Kosterlitz-Thouless |33| vortex unbinding transition. Even below T c 



however, there is no privileged direction in the tilt field since a broken continuous symmetry 
is impossible in two dimensions for systems with short range interactions. Thus, one 
expects the polymer to crumple for T < T c , and one can ask for the wandering exponent 
v governing the mean end to end distance, R ~ N u , where N is the number of monomers. 
Note that v = | in the three-dimensional problem discussed above. The problem for d = 2 
is nontrivial because below T c , the nematic field has long range correlations decaying with 
a continuously varying exponent: 

(n(x) ■ n(x')> = (e^We^M) ~ |x - x'P (T ) . (2.9) 



The exponent r](T) varies from ?] = OatT = Oto?7=4;at the transition |3^]. There are a 



priori two Frank constants K\ and K% in two-dimensions, but it has been shown j34] that 
at large length scales two dimensional nematics become isotropic (K\ = K3 = K) and 
are described at long wavelengths by the usual XY model free energy § fdx (V(j)) 2 . For 



T < T c one has r\ = T](T) = Above T c the correlations decay exponentially in Q2.9|) , 

and then one expects that the wandering exponent of the polymer is the pure self-avoiding 
random walk value v = 3/4 (or v = 1/2 for an ideal chain) . 

A simple random walk argument gives an interesting prediction for the exponent v 
below T c . In the infinite g limit one expects that the polymer be totally aligned with the 
local field. 

^ - e l ^ s » , (2.10) 

where we have used the complex notation for the position z(s) = xi(s) + 2x2(5), s being 
the arclength along the polymer. Upon integrating over s we find 

(\z(s)-z(0)\ 2 ) z{0U = f f\udu\e^ z ^-^ u '^ z{0U (2.11) 
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where the average is over both initial conditions z(0) for the polymers and over the 
smoothly varying director angle 4>(z). We now neglect the reaction of the polymer on 
the nematic field and replace the correlation in the right hand side of ( |2.11| ) by the cor- 



relation of the unperturbed nematic field (|2.9|) . This procedure amounts to replacing the 
problem by the one of determining the wandering of the tangent curves of a pure nematic 
(see Figure 4). We are led to an integral equation, namely 

(\z(s)-z(0)f) z{0) = f fdudu't C ° nS ) \ (2.12) 

Jo Jo \l¥)-¥')r/ < (o) 

Upon assuming that \z(s) — z(0)\ scales as s u , we find that the self consistent value of v is 

u(T) = 2 —— . (2.13) 

Our results coincide with those of Flory theory applied to random walks in quenched 
random environments, which has been argued to be exact for divergenceless flows [|3~2"| , |3"5|| . 

According to Eq. ( |2.13j ) v is always close to v = 1, continuously decreasing with 
increasing temperature from v(T = 0) = 1 to u(T = T c ) = 8/9. Note that v always 
exceeds 3/4, the value of a self-avoiding random walk in two dimensions. Our neglect of 
self avoidance for T < T c is thus self-consistent because of the relatively small number of 
self-intersections which occur for v > 3/4. 
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2. 3. Field Theoretic Treatment of Hairpins 



In this section we study a single polymer embedded in a rigid nematic matrix, and 
show that hairpins induce an Ising-like crumpled state. A field theory due to Cardy ]2^] is 
used to describe the polymer and introduce hairpins, providing a simplified example of the 
boson field theory used in later sections. We show explicitly that hairpins cause a single 
directed polymer to crumple into a more isotropic configuration at long length scales, as 
suggested in |J . Although crumpling always occurs in a single polymer, we show in Section 
7 that interactions induce an Ising-like transition from a dilute crumpled state to directed 
state stabilized by interactions when many polymers are present. 

Consider the following model free energy for a single polymer in a nematic matrix n, 



ds 




9 I r ' 
- [T ■ n 

2 



(2.14) 



where T 



dR(s) 



is the three dimensional tangent vector to the polymer, and s is the 
arclength. The first term is the bending energy, while the second represents the lowest 
order coupling to the background nematic compatible with overall rotational symmetry. 
Higher order terms in T are possible, but are inessential for the following discussion. Note 
that \f\ = 1. 

Suppose the fluctuations in n are suppressed by imposing a large magnetic field, so 
that n = z. Then one has 



ds 
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(2.15) 



which is exactly a one dimensional Heisenberg model with a quadratic Ising-like anisotropy. 
As a function of T z , there are two minima for T z = ±1. If g is very large, one can expand 
perturbatively around each minimum. Because this is a one dimensional problem, the 
symmetry is in fact restored by tunneling events between the two minima, which correspond 
to hairpins. T z thus plays the role of an Ising variable and hairpins are analogous to an 
Ising domain wall along the one dimensional chain. The tunneling probability from one 
minimum to the next is the energy of a hairpin occurring over a distance R, with energy 
eh ~ minR + ~ (g^) 1 / 2 , and the radius of the hairpin is Rh ~ ( K /fl r ) 1 ^ 2 which 
coincides with de Gennes more elaborate calculation || . Note that the size of a hairpin is 



also of the order of the Odijk deflection length A [10]. The typical distance between two 
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hairpins along the chain is thus / ~ ae th l k ^ T 1 which can be very large. The statistics of the 
hairpins is analogous to the statistics of one dimensional Ising domain walls (with Ising 
coupling J ~ eft,). This analogy implies that the length scale a is just the Odijk deflection 
length. Note that a stretching force applied to the endpoint of the chain is equivalent to 
an Ising magnetic field, since j ds hT z = h(z(L) — z(0)). The Ising magnetization is the 
analog of the total size of the polymer along z. The elastic modulus G for Hooke's law 
along z is thus given by the Ising susceptibility, G ~ (l/T)e J//fcBT as T — > 0. 

The effect of director fluctuations will be discussed in Section 7. Here, we continue 
to neglect them and implement the above ideas with an Ising-like field theory. As noted 
by Cardy the propagator of a single directed polymer of variable length along z with 2 
transverse dimensions can be written as the correlation function of a quadratic action, 
which we write directly in the continuum limit [^9| , 

G(r, 0;z, 0) = i J ViJjViP* ip*(0, 0)^(r, z)e~ So (2.16) 

where 

So = J dzd 2 riP* (d z - DV\ - p) tp (2.17) 

with D = k B T/2g. Here G(r, 0; z, 0) = E W N (r, 0; z, Oje - ^ and W n {t, 0; z, 0) is the total 
weight for directed walk of iV steps to begin at (0,0) and end at (r, z). The chemical 
potential ft must be adjusted to give the correct polymer size. The use of a retarded prop- 
agator only and the neglect of self-interactions if hairpins are excluded has been justified 
by Cardy. The critical point of this quadratic theory is p, — > 0~ corresponding to infinite 
length chains. The propagator in Fourier space is then 

G(<U, q z ) = . * 9 = G (^, q z ). (2.18) 

Now we discuss a single directed polymer with hairpins but without self avoidance 
(see Figure 5). One can introduce hairpins by simply adding to the action (|2.17|) the term 



s = s ° + f / dzd2r V 2 + (2 ' 19) 



which allows for "pair creation". Let us consider the new term as a perturbation and 
expand in w the above correlation function ( |2.16| ). One then generates diagrams as in 
Figure 5, with a factor w per hairpin and Go(q 1 ,g z ) per solid line. All closed loops are 
canceled by the normalization factor Z, so the n —>■ trick is unnecessary. The term of 



11 



order w 2n is the sum of all possible ways to go from (0, 0) to (r, z) with 2n hairpins - the 
odd terms vanish. The analogy of hairpins with kinks in an Ising-like system shows that 
we should take w oc e £h / fcBT , with eh ~ y/g^- If one follows the connected line from (0, 0) 
the factor associated with the part which represents backward propagation (e.g., after an 
odd number of hairpins) is actually Gq((\ ± , q z ). Thus the sum of all these diagrams is 

[l+(|) 2 Go^ + (f) 4 (G ^) 2 + .. 



G = G ((i ± ,q z 
Gq 



(2.20) 



1 ^(-TOCtq 

Thus the propagator is now 



G= !^±M^g _ (2 . 21) 
Q 2 z + (Dqi-fi) 2 -f 

The critical point has been shifted and is now at Jx = —w/2. This corresponds physically to 
the fact that allowing hairpins adds an extra entropy in the system, and the proliferation 
of paths occurs earlier when \x is raised from — oo. To find the asymptotic large distance 
behavior, we set ft = —w/2 + dfi and expand for small 5 ft, q^, and q z . Then we have 

° ~ (<fc/«0 2 +Dql ~ Sfl (2 ' 22) 
Note that complete propagator is actually twice the above result (more precisely it is 
G + G*) because one has now to allow for graphs where the first propagator in the above 
series goes backwards. This is the propagator for an (anisotropic) Gaussian random walk 
in (i-dimensions. The polymer size now scales as VW both along z and along r, but with 
an anisotropic radius of gyration tensor. The directed nature of the walk has disappeared. 



3. Field Theory for a Liquid of Chains with a Nematic Background 

In this section we use the analogy with the statistical mechanics of two-dimensional 
bosons to describe the properties of chains. The details of this analogy are discussed in 
|l],|36|], so we will indicate here only the differences with the case of the flux lines. For a 
fixed configuration of the nematic solvent, the analogy works as before, except that now 
the bosons are also interacting with an external "time" and space dependent field. The 
average over configurations of the nematic field has to be carried out at the end. We 
initially discuss the theory for chains which span the system. Internal free ends are then 
introduced by adding a source to the boson coherent state field theory. The field theory 
with the free ends is then used to calculate some important correlation functions in Section 
4. 
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3.1. Mapping onto two-dimensional bosons 

We start from the following partition function for the N chains in a nematic field 
discussed in the introduction 

N 



i=l 



Z N [Sn} = ± J f[Vr^ 

N 



i=i 



r iz , z 



Y + f2V(r t (z)-r J (z)) 



(3.1) 



where we have omitted the bending term and restricted our attention to interactions be- 
tween lines through an equal "time" potential V{ri{z) — Tj(z)). The quantity which we 
are ultimately interested in is the average: 



V5n(r, z) Zn[5ji] exp 



KT 



(3.2) 



where F n is given by 

Well known transformations [|37| , |3~8|| using transfer matrices allow to write Zjv[5n] in 
the Hamiltonian form (first quantization) , 



Z N [5n] = J rf2r i • • • d 2 r N d 2 r[ . . . d 2 r' N 



r'... r' 



N 



Texp 



dz 



n{z) 



ri . . . rjv 



(3.3) 



where |ri . . . rjv) and |r^ . . . r^) are states corresponding to the entry and exit points of the 
polymers |I[ and the time ordering operator T is necessary a priori since Tt is z dependent. 
H can be deduced from the Lagrangian 



i i<j 



(3.4) 



according to the usual rules of Euclidean quantum mechanics [|3"?| 

H(p k , rjb, z) = C + i ^ pj ■ Tj , (3.5) 



where Yj has been eliminated using pj = i^-. The quantization rule is then pj = 
—ikvTVj. We define V so that it operates only within the plane perpendicular to z. 
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The quantity r will always mean a vector perpendicular to z. Following these rules, we 
obtain the Hamiltonian for a fixed configuration of 5n: 



H = J2 V ? + k - T E \(yM^ *) + *n( rij *)Vi) + ^(r, - p,-) , (3.6) 



where we take a symmetrical ordering []37j . 

A coherent state representation can now be developed in analogy to the treatment of 
flux lines in Pj] , and one obtains for the grand canonical partition function: 

oo „ 

z) exp — S[ij)*,ij), 8n] , (3.7) 

N=0 J 

where ijj(r, z) is a complex boson field. The boson action S reads 



S 



dz l d\ 



^{r,z) (J--DVl-jij 1>(t,z) 
+ ^ [ip*(r, z)V_L^(r, z) - ip(r, ^)V±^*(r, z))} ■ 5n(r, z) 



+ 



+ - J dr' v (r — r') \ip(r, z)\ 2 \ip(r' , z) 
F n [Sn] 



(3. 



Similar manipulations show that the density of flux lines is 



p(r,z) = |V>(r, z)\' 



(3.9) 



where D = k B T/2g, ft = p/k B T and v = V/k B T. In Section 4, we shall calculate corre- 
lations in the density p(r, z) = 52 i=1 8[r — ri(z)] may be calculated via the identification 
( p.9|) . In the mean field approximation, we have po = (\i(j\ 2 ) = ft/ v. The field theory 
embodied in ( |3.8|) differs from the action for flux lines |l|] only in the coupling of the boson 
"current" tjj*V ±ifj — ±ifj* to the director field. As discussed in Appendix A, rotational 
invariance of the original "Lagrangian" model (|3.1| ) (i.e. "Lorentz invariance" of the fic- 
titious bosons) forces the coefficient of (■0*V X '0 — i^V ± ip*) ■ Sn to be exactly half that of 
i\)*d z i\) in the second quantized coherent state representation. 
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3. 2. Effects of Finite Chain Length 

To allow for polymers of finite length, which start and stop in the interior of the 
sample, we add a source term to ( ft.8| ) 



S -> S -h J dz J ' cfrty + ip*]. (3.10) 



Upon replacing ip and tp* by their mean field average value y/po, we see that the quantity 
hy/p~o is the probability per unit area and per unit "time" z of starting or terminating 
a polymer. Since an average of one polymer will thread a cross-sectional area p^ 1 per- 
pendicular to z, the typical polymer length associated with this Poisson-like process is 

A more formal proof of the equivalence of bosons with a source to the statistical 
mechanics of directed polymers of finite length can be constructed along the lines taken 
for isotropic polymer melts |39|,f|0] . We first expand the partition function associated with 
( j3.8|) in h, v, and the nonlinear coupling to the director field. The Fourier transformed 
propagator in the resulting Feynman diagrams (see Fig. 5) is 



or in real space 



G (qx,q z ) = — - 2 , (3.11) 

-iq z + Dqj_ 



G (r, z) = (^(r, z)r(0, 0)) = S^-^—e^^*, (3.12) 

4nDz 



where Q(z) is the step function. The expectation value in ( |3.12|) is taken with respect to 
the Gaussian action 

S = J dz J d\ [ip*(d z - DV 2 ± )tP] (3.13) 

Equation (|3.12|) is identical to the propagator used in the polymeric description of flux 
lines presented in 0. We can think of ^(O, 0) as creating a polymer at (0, 0) and i/j(r, z) 
as destroying this polymer at (r, z). Although this produces the usual diffusive random 
walk propagator for z > 0, Gq(t, z) = when z < 0, showing that propagation backwards 
in the time-like variable z is impossible. 

The grand canonical partition function may now be expressed formally as 

oo oo 

Z ^= E ^ZfrjD^K^h* (3.14) 

N m =0p=0 
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where Z V N (D,v,Ki) is the partition function for N m monomers distributed among p 
polymers and L m oc N m is the total length along z occupied by the N m monomers. Note 
that h 2 plays the role of a polymer fugacity, while e M controls the density of monomers. 
Figure 6 shows a typical contribution to Z V N of order h 8 . The solid lines are the polymer 
propagator Go(q ± ,q z ) discussed above. Dashed lines represent the interaction potential 
V(r), while the dotted lines signify interactions induced by the background nematic field. 
The retarded nature of the propagators insures that large numbers of "unphysical" graphs 
disappear. The graph shown in Figure 7a, for example, with zero momenta on its external 
legs, is proportional to 

dq z d 2 q ± |U(qJ| 2 
(2tt)3 (- iqz + Dq± 2)2> \ ■ > 

which vanishes because both poles in the q z integration are on the same side of the real 
axis. The graph shown in Figure 7b does not contribute for similar reasons. The graphs 
shown in Figure 7c are constants which can be absorbed into a redefinition of the chemical 
potential. 

By using the expansion ( |3.14| ), we can easily show that the average polymer length is 
given by 

(L m ) = ^lnZ gT = ^ } (3.16) 

while the average number of polymers is 

(p> = „_!_ lnZgr = E^ES^JA^W* {317) 

We can now calculate the typical polymer length 

i = {L m )/(p) 

= (Mfezjf) (3-18) 
|(^(r,z)+V*(r,z)) 

directly from the mean field approximation to the partition function ( |3.7j ). We assume for 
simplicity the contact potential V(r) = k B Tvd 2 (r). Upon making the substitution ( |3.10| ) 
in ( p.7| ) and setting ip(r, z) to a constant value ip = = = J~p$, we have 



lnZ gr = — fimin | —[i^ 2 + -vi/jq — 2hipQ \ (3.19) 

where O is the (three dimensional) volume. There are two limiting cases to consider. When 
JJl ^> 0, we assume an ordered state only slightly perturbed by the small source field h. 
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As we shall see, this means restricting our attention to very long chains. The minimum of 
( p,19| ) is then given to lowest order in h by 



V>0 



1 + 



v h 




fi 3 2 



(3.20) 



while the partition function is 



Z gr = exp 



2v V v 



(3.21) 



We then find from ([PJD and fl3l7D that 



and 



so that a typical length is 



(L m ) « ^ 



(p) « Vo^ 



(3.22) 



(3.23) 



(3.24) 



as also follows directly from the mean field approximation. Note that this length diverges 
as h — > 0, and agrees with the estimate made at the beginning of this section. 

If /i<0, the polymers are dilute, we are in the disordered phase of the model. We 
assume that polymers follow an imposed external direction, as in electrorheological fluids, 
so that we need not worry about the crumpling transition discussed in the Introduction. 
One can then neglect v and find that the order parameter is 



tpo ~ h/\n\ 



(3.25) 



while the grand partition function is 

Z gr w exp 
The total length of polymer is now 



(3.26) 



while the average number is 



h 2 

\L m ) = 



h 2 

(P) = ^ 



(3.27) 



(3.28) 
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A typical polymer size is thus 



1 

lAl 



(h = o+) 



(3.29) 



which increase as the transition is approached from negative p. 

Experiments are usually done varying the polymer concentration with a fixed distri- 
bution of polymer sizes, and hence, fixed I. This produces the trajectory shown in Figure 8 
on the phase diagram plotted as a function of —jj, and ifto. Typical polymer configurations 
at two points on this phase diagram are indicated. Note that the model corresponds to 
a distribution of chain lengths instead of a monodisperse sample. This is also a feature 
of the des Cloiseaux model of isotropic polymer melts, for which it is possible to draw a 



very similar phase diagram |39|,(i0|]. A broad distribution of chain lengths is probably a 
good approximation for electrorheological fluids and ferrofluids, in which the chains are 
constantly breaking and reforming. The difference between polydisperse and monodisperse 
samples is not, in any event, expected to be important when the average chain length is 
large ||0 . 



4. Results for Correlation Functions 

l — 1/2 

We now assume that the polymers are dense and entangled, i.e. yDl ^> p 1 . 
Correlations for small h can then be calculated as in fl],|36| , by expanding about the mean 
field order parameter ijjQ = ^/po « y/p,/v[l + O(h)]. 

4-1. Gaussian Form for the Action 

Let us write the complex field if)(r,z) in terms of the density p(r, z) and the phase 
6>(r, z) as: 

-0(r, z) = p(r, z) 1/2 exp(z#(r, z)) . (4.1) 

The measure takes the simple form, up to unimportant constant factors: 

£V*(r, z)Vip(r, z) = Vp(r, z)V9{r, z) . (4.2) 

To simplify the discussion, let us again consider a contact potential V(r) = Vo<5 2 (r), and 
set v = Vo/k B T. The action then takes the form in the new variables: 

D 



S = J dzd 2 r 



i P d z e + —(v ±P ) 2 + Dpiwj) 2 + i P v ± e ■ Sn 

Ap 

ftp + -vp 2 - 2hp 1/2 cos 9 
2 



(4.3) 
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Upon setting p(r, z) = po + Sp(r, z), where po = ifi 2 is given by Q3.2Q ), and expanding to 
quadratic order the fluctuations dp and 9, the action becomes, up to total derivatives and 
constants, 



S = dzd r 



iSpd z 9 + ^-{VJpf + Dp (VJ) 2 + i Po V ± e ■ 5n 



(4.4) 



Note that the term involving 5n can be rewritten as —ipo9(V ± ■ <5n), so that only the 
longitudinal part of Sn couples to the polymer degrees of freedom. 
Upon expanding the fields in Fourier modes 



9{v,z) = (O)" 1 / 2 e lzq * +tcl± - r 9(ci ± ,q z 



(4.5) 



q±,qz 



with similar expansions for 5p and 8n, one obtains the action: 

S= \ X+{^q z )G- l {^q z )X{^q z 



(4.6) 



where X = (9, 8p, Sn x , Sn y ) and, (taking H = in ( |1.1| ) for simplicity), 
G~ l = 



2Dp q ± 
+2hp\ /2 



-POQx 



-poq y 



2 P0 



POQx 



POQy 



+hp Q 






3/2 



(K iq 2 +K 2q 2 y +K 3 q 2 ) Kl _ K2 



k B T q * q y 



k B T qxq y 

(K 1 q 2 ! +K 2q 2 +K 3 q 2 ) 

fc R T 



(4.7) 



It is now straightforward to calculate any desired correlation function, by inverting this 
4x4 matrix. 

4-2. Discussion of Correlations 
The structure function 



S(<l ± ,q z ) = (\5p(c ix ,q z )\ 2 ) 
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may be written 

S(a a) = n m l + ^(q) (T 1 + Dgj) 

iqx,(/2j W ^(g)[g z 2 + e2( ? J] + ipog!6-2( ?x )/(£-i + J Dgi ) ^ • ^ 

where 

^ = ^I±^ (4.10) 

we have set p\j 2 lh = £, and 

e 2 ( ?x ) = (r 1 + ^^(r 1 + £>gi + 2p oW ). (4.11) 

To obtain results for arbitrary in-plane pair potentials V(r), let v — > U(q x )//c B T, where 
t^(qj_) is the Fourier transform of U(r). Unlike its flux line counterpart, the "Bogoliubov 
spectrum" ( |4.11[ ), has a gap as q ± — > 0, due to the finite polymer lengths. Upon evaluating 
(Sni(q)5rij(—q)), we determine the renormalized Frank constant via the identification 

(SmffiSnji-q)) = -—L—fPfi + 1 2 Ph (4.12) 

K<i and K% are unrenormalized, while 

*i*(<U, ».) = + ^r gg^—L-^ . (4.13) 

Upon taking the limit K(q) —>■ oo, which suppresses all fluctuations in the director 
field, ( |4. y| ) can be applied to polymer melts with an externally imposed average direction. 
In the limit of small wavevectors, the structure function takes the simple form, 



pjql + K/G 
Bq\ + Kq 2 z + G~ x KBp^ 



S(^q z ) = k B T — - ,To - 2 (4-14) 



where B = PqVq + 0{£~ x ) is a two dimensional bulk modulus, K = p g is a tilt modulus 
and G is related to the average length via G = £k B T/2p Q . When £ — > oo the structure 
function contours are straight lines passing through the origin, with no scattering along 
the line q ± = |||7| . Scattering reappears at q ± — for finite £, however. Fits of (|4.14| ) 
to scattering data should allow a simple direct determination of the average chain length 
£ and the important parameters B and K. 

The small wave vector limit of ( ^.9| ) when the Frank constants are finite has a different 

form 

S(a a)-kT plql + {Kiql + K * ql)/G (4 15) 

6(<U,fcJ - k B T Bq[ + {BjGpl + + R ^ . (4.15) 
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When G = £k B T/2p — > oo, QH5Q reduces to the prediction |31[] of a hydrodynamic theory 



due to de Gennes JJ(J. Finite length effects distort the contours near the origin, however, 
and are likely to be quite important in fitting real scattering date. Fits to ( f4.15| ) should 
lead to direct experimental determination of the Frank constant, the bulk modulus B and 
the mean polymer length. Upon taking the limit q — * in we find the dependence 

on the renormalized splay elastic constant -fTf on the polymer length, 

K? = K X + h B T£ Po , (4.16) 

in agreement with a prediction of R.B. Meyer ||, but in disagreement with a suggestion 
by de Gennes ||. 



5. Hydrodynamic Treatment of Correlations 

The hydrodynamic description of isotropic liquids of atoms or small molecules has 
been understood for many years |[41|| . The long wavelength density fluctuations are Gaus- 
sian, and static correlation functions are described by Ornstein-Zernike theory |4!|. We 
generalize here an analogous theory of liquids of oriented lines in three dimensions. The 
basic concepts were first discussed by P.G. de Gennes and R.B. Meyer in the context of 
polymer nematic liquid crystals over a decade ago []8]|J. These ideas were recently used 
to determine the form of the liquid crystal structure function near the origin of reciprocal 
space by Selinger and Bruinsma [BT|. The hydrodynamic theory can in fact be derived di- 
rectly from the boson representation, as shown in Appendix B. We show here how a finite 
density of chain lengths can be incorporated in a natural way for ferro- and electrorhe- 
ological fluids, as well as for polymer nematics. The results agrees in all cases with the 
long wavelength limit of those in Section 4. We also use the hydrodynamic theory to show 
that there are essentially no differences at long wavelengths between polymers in a nematic 
solvent and dense polymer nematics in an isotropic solvent. Some of our conclusions were 
reviewed recently in [[13] . 



Hydrodynamics will, in addition, allow us to discuss the physical meaning of the boson 
order parameter introduced in Section 3. To this end, we calculate the energy of an isolated 
free end in a ferro- or electrorheological fluid. Such a calculation has already been carried 
out for polymer nematics by Selinger and Bruinsma [|3 . 
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5.1. Ferro- and Electrorheological Fluids 

We first assume for simplicity that the chains of magnetic or electric dipole particles 
span the system along the z-axis. The effects of finite chain length will be put in later. 
The basic hydro dynamic fields are now an areal particle density 

N 

Pm ic(r,z)=Y,S 2 [r-r J (z)] (5.1) 

3=1 

and a "tangent" field in the plane perpendicular to 5, 

W c ( ri z)^^ 2 [r-r 3 ( Z )]. (5.2) 
i=i 

We coarse grain these microscopic fields to obtain smoothed density and tangent field 
p(r, z) and t(r, z). 

We now expand the free energy of the liquid to quadratic order in the density deviation 
Sp(r, z) = p(r, z) — po and in t(r, z) 

F = — j / d\dz [B(5pf + K\t\ 2 ] . (5.3) 

The parameter B is a bulk modulus for areal compressions and dilations perpendicular to 
the z-axis, while K is the modulus for tilting lines away from the direction of the applied 
field. Because we are dealing with lines, and not simply oriented anisotropic particles, 
( p73|) must be supplemented with an "equation of continuity," 

d z 8p + V_L-t = 0, (5.4) 

which reflects the fact that vortex lines cannot stop or start inside the medium. Corre- 
lation functions can be calculated by assuming that the probability of a particular line 
configuration is proportional to exp(— F/k B T), and imposing the constraint ( |5.4| ) on the 
statistical mechanics. 

We now modify ( |5.3| ) to allow for chains which stop and start inside the medium. If 
the chains are long and entangled, we can treat the chain heads and tails as independent 
ideal gases, following similar ideas for polymer nematics by R.B. Meyer 0. The free energy 
(|5.3|) is now replaced by an expansion of the form 



F = - I ifrcli 



2 2 

B[¥) +k(—) +G(d z Sp + V ± -t) 2 
PoJ \PoJ 



(5.5) 
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which differs from only in a term proportional to the square of the constraint displayed 
in ( |5.4|) . Although other terms proportional to gradients of 5p and t can appear, we have 
kept only those couplings which dominate in the limit of long chains. We can now treat 
Sp and t as independent fields and then take the limit G —>■ oo to impose the constraint. 
The coupling G is finite, however, when the chains are of finite length. To determine the 
value of G in this case, note that the chain heads and tails act as sources and sinks on the 
right hand side of the conservation law (|5.4|) . It follows that || 

d z 5p + V_L-t = p H - p T , (5.6) 

when finite densities of chain heads p H (r,z) and chain tails p T (r,z) are present. If the 
heads and tails are treated as two noninteracting ideal gases, a term of the form 

6F = -G [ d 2 rdz (p H - p T f (5.7) 



2 

should appear in the free energy. The coefficient is just the concentration susceptibility 



for an ideal binary mixture, which is well known to be [[45 



G = I f f / \ ( 5 - 8 ) 
\Ph) + {Pt) 

where (p H ) and (p T ) are the average concentrations of heads and tails, respectively. Now, 
each chain contributes both a head and a tail, (p H ) = (p T ) = p c hain, and the three- 
dimensional density of chains is p c hain = Po/^ where I is a typical chain length. It follows 
that 

G=^ (5.9) 
2po 

which diverges as I — > oo. It is easy to check that the structure function which results 



from this hydrodynamic treatment agrees with the result (|4.14|) of the more microscopic 
boson calculations of Section 3 in the long wavelength limit. 

Figure 9 shows the scattering contours expected for dense aligned ferro- or electrorhe- 
ological fluids. The maxima along the </ x -axis occur approximately at the position of the 
first reciprocal lattice vector in the nearby crystalline phase, and are not accounted for by 
the hydrodynamic theory. As discussed in [|l|, the half- width at half-maximum along q z 
for fixed q ± controls the decay of density fluctuations along the z-axis. Hydrodynamics, 
(in agreement with the boson theory) determines via ( f4.14j ) the linear form of the contours 
near the origin. The rounding of these contours (indicated by the dashed lines) due to 
the finite polymer length is one of the principal predictions of [2(| and this paper. When 
£ — > oo, the contours remain linear as g ± ,g z ^0 and the scattering vanishes along the 
g z -axis, as in the case of flux lines ffl . 



23 



5.2. Polymer Nematics 

The hydrodynamic treatment of infinitely long polymer nematics is due originally to 
de Gennes ||30|| , and was applied to correlation functions by Selinger and Bruinsma [|31| . We 



first assume an isotropic solvent. The hydrodynamic variables are now the areal polymer 
density ( |5.1[ ) and a fluctuating nematic deviation field <5n(r, z) attached to the polymers. 
We again allow for a dilute concentration of chain ends by writing the hydrodynamic free 
energy as 



d\dz 



B [ — ) +G(d z 5p + Po V ± -5n" 2 



+ F n [5n] (5.10) 



.Po. 

where F n is given by (|1.1|) , and G must again be given by ( |5.9|) in the dilute limit. When 
G —>■ oo, we recover the constraint, 

d z 5p + p X7 ± -5n = (5.11) 

For finite G, one again finds agreement with correlations calculated from ( |5.10| ) and the 



hydrodynamic limit (|4. 15|) of the result ( [Op obtained from the microscopic boson theory. 



Figure 10 shows the hydrodynamic predictions for scattering off the polymers in the 

— 1/2 

limit of small momentum transfer. The structure for larger q ± ~ p should be similar 



to that shown near the peaks in Figure 9. ||25|| . For infinitely long polymers, the scattering 
vanishes along the g z -axis, and the contours take the form q z oc q ± Pl[] . The boson and 
hydrodynamic theory prediction (|4.15| ) leads to the rounding indicated by the dashed lines 
in Figure 10, an effect which is likely to be quite important in fitting real experimental 
data. 



As discussed in P3(j , hydrodynamics also makes interesting predictions about freeze 
fracture experiments on directed polymer melts. It can be shown, in particular, that 
density fluctuations measured in a fracture plane perpendicular to z provide a precise 
signature that one is dealing with a liquid of lines rather than a liquid of points. Under 
favorable circumstances, it is even possible to determine a typical polymer length from 
such measurements. 

We can also treat polymers in a nematic solvent using hydrodynamics, and check 
that there are no significant differences with conventional polymer nematics in the long 
wavelength limit. The (unpolymerized) nematic solvent will now be described by the free 
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energy (LI) and we introduce coarse grained polymer variables as in (|5.1|) and ( |5.2| ). The 
hydro dynamic free energy (as derived explicitly in Appendix B) is then 

2 



B 



dp 
Po 



+ g\t - poH 2 + G(d z Sp + V ± -t)< 



+ F n [5n] 



(5.12) 



where the coupling proportional to g reflects the analogous coupling in (|2.1|). We now 
integrate out the solvent degrees of freedom, which leads to an effective free energy 



F e ' ff /fc B T = / V s ne -F'/k B T 



given by 



where 



Ft 



eff 



d 2 rdz 



^ ) +G(8 z 5p + V ± -t) 
Po. 



+ K[t] 



(5.13) 



(5.14) 



KM 



d d q ± dq z 
(2n) d 2tt 



Kiqj + K 3 qj gjgj 
K iq 2 ± + K 3 q 2 z + gpl q\ 

K 2 q 2 ± + K 3 q 2 z & - q^qi 



9 



+ 9 



ti(<f)tj(-q) 



(5.15) 



K 2 q 2 ± + K 3 q 2 z + gpl q 2 ± 
At long wavelengths, the coefficients of the longitudinal and transverse projectors in ( |5.15| ) 
simplify, and F' n for polymers in a nematic solvent reduces to the F n appropriate for 
conventional polymer nematics. The two systems are indeed equivalent from the point of 
view of long wavelength polymer correlation functions. 

5.3. Defect Energies and the Boson Order Parameter 

Hydrodynamics also allows us to better understand the boson order parameter used 
in Sections 3 and 4. For an analogous discussion for flux lines in high T c superconductors, 
see 0. Note that the representation (4TT) only makes sense if there is "phase coherence" 
in the equivalent "boson" system. Consider the correlation function 



G(r,r';z,z') = (i;(r,z)r(r',z')} 



(5.16) 



where i/j(r,z) and ?/>*(r', z') are, respectively, creation operators for polymer heads and 
tails. We assume that the polymers are long and entangled, so that h = in ( |3.10| ). Phase 
coherence means long range order in G(r, r'; z 1 z'), 



lim 

|r— r' I— >oo 



lim G(r, r'; z, z) = const > 0. 



(5.17) 
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To understand what this long range order means, consider first the behavior of ( p.!6|) in a 
hexagonal crystal of directed polymers, with the strands again aligned with the z-axis. The 
composite operator in ( |5.16|) creates an extra line at (V, z'), (i.e., a column of interstitials 
in the solid), and destroys an existing line at (r, z), creating a column of vacancies. As 
shown in Figure 11, the lowest energy configuration is then a line of vacancies (for z' > z) 
or interstitials (for z' < z) connecting the two points with an energy as proportional to 
the length s of this "string." The string tension a will depend on the angle 9 this line 
makes with the z-axis. It follows that the correlation function ( |5.16|) decays exponentially 
to zero (i.e., like exp(— a(9)s/k B T)) for large separations in this crystalline phase. In a 
directed polymer melt, on the other hand, the concept of vacancy and interstitial lines has 
no meaning, and the string tension a will vanish for large s, implying long range order in 

G(r,r';z,0- 

Long range order in the boson order parameter ijj(r, z) means that 

(V>(r, z)) = (V>*(r, z)) oc e ~ E ' /kBT > 0. (5.18) 

Here, E* is the energy of an isolated polymer head or tail. The hydrodynamic theory pro- 
vides a transparent demonstration that unlike in hexagonal crystals, this energy is finite. 
This energy has already been calculated by Selinger and Bruinsma for polymer nematics 



4JJ], so we concentrate here on the hydrodynamics for ferro- and electrorheological fluids. 
The constraint ( |5.4|) applies everywhere away from an isolated head or tail, and is conve- 
niently implemented by expressing Sp and t in terms of a two component "displacement 



field" u(r,z) 



dp = —poV ± ■ u (5.19a) 

t = Po^ (5.196) 

oz 

following a similar trick by Taratura and Meyer ||6| for polymer nematics. The free energy 
becomes 

F = ^ J dzd 2 r [B(V ± ■ u) 2 + K(d z u) 2 ] . (5.20) 

The extremal equations associated with ( p.20| ) for an isolated polymer head at the 
origin are 

SV x (Vx ■ u) + Kd 2 z u = B P q 1 6(z)V i o 2 (r), (5.21) 
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where po is the average in plane polymer density and @(z) is the step function, Q(z) = 
1,2 > and O(z) = 0, z < 0. See (47] for a discussion of the source term for the closely 
related problem of a dislocation in a two-dimensional smectic liquid crystal. Here, the 
source represents the absence of a line along the positive z-axis. In a crystal, this line 
would be a string of vacancies, and its energy would be infinite, due to the disruption of 
the crystalline order parameter in the vicinity of the line. In a liquid, however, the vacancy 
free energy vanishes and we need only consider the long range strain field associated with 
( |5.21|) . Upon solving ( |5.21|) in Fourier space, we find 

q z {Bqi +Kqj) 

from which leads to an explicit expression for the tilt field 

B 1 / 2 K r , . 

t = Pod.u = 4% {Kr2 + Bz2)3/2 - (5-23) 

The density perturbation associated with the free end is 

ip =- P0 (V i .u) = Z^_i_. (b.24) 

Since both these "strains" fall off like l/(distance) 2 far from the origin, the integrated 
strain energy in ( |5.18|) is indeed finite, i.e., 

E* < oo, (5.25) 

where we have included a microscopic short range contribution in the defect energy E*. 
The analysis for polymer nematics [11] is more complicated but also leads to the con- 



clusion that the energy of an isolated head or tail is finite, consistent with a nonzero value 
of the boson order parameter. The strains in ( |5.23| ) and (|5.24p resemble those expected 
for a magnetic monopole in an anisotropic medium. Indeed, ( [5.4] ) is just the condition of 
"no magnetic monopoles," V • b = 0,if we identify Sp with b z and t with b^, while (|5.3|) is 
just the magnetic field energy of a medium with an anisotropic permeability. The energy 
of isolated magnetic monopoles is indeed expected to be finite in three dimensions, except 
when the field lines form an Abrikosov flux lattice [EL 
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6. Renormalization Group Analysis of the Dilute Limit 

The analysis in the previous sections applies when the chains are dense and entangled. 
To treat the dilute limit, we perturb about the limit in which (if)) « 0. To this end, we again 
consider the (coherent state) functional integral representation of the partition function: 

Z gr = JvipThl>*V6nexp-S[ip,ip*,6n], (6.1) 
and break the action into three parts, 

S[iP, V>*, Sn] = J d d r J dz [ijj* (d z - DV\_ -p)^ + t#| 4 ] 



+ I d d r I dz8n-(ip*V±ip-ipV±ip*) 



(6-2) 
+ F n [6n]/k B T 

This form of the action is just ( |3.8j ), specialized to the case v(r) = vS d (r) and with a new 
parameter A to help organize perturbation theory in the coupling between the polymer 
and nematic degrees of freedom. As in in |l|, we impose a cutoff A on the perpendicular 
wavevectors q^, of order the inverse range of the interaction. For now we neglect the 
possibility of free ends and do not introduce a source term. It is convenient to consider d+1 
dimensional directed polymers with <i-directions perpendicular to the average direction. We 
take the limit d = 2 at the end of the calculation. 

We can now attempt to expand quantities of physical interest in the nonlinear cou- 
plings v and A. Consider, for example, the propagator, 

G(r,z) = (V(r,^*(0,0)), (6.3) 



or, in Fourier space 



This can be written as 



1 



(6.4) 



(6.5) 



-iq z + Dql - fi + S(q_L, q z ) 
where, to lowest, non-trivial order in v and A, the self energy graphs shown in Figure 12 
lead to 



£(q^,<? 2 



dk 7 fd d k 



(2<L 



2tt 
- k 



{27T) d 

(2^- 



- 2v- 



k z ) + D(q x -k 1 _) 2 
1 

-ik z + Dk\ — \x 



8_v - k\k\jk\ 
fx \ K x kl + K 3 k 2 z ' K 2 k 2 ± + K 3 k 2 z 



+ 



(6.6) 
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The first term comes from the nematic-polymer interaction, while the second term arises 
from the polymer-polymer interaction. This self interaction of a single polymer is not 
present in the original microscopic model ( |1.2| ), and arises here because of the nonuniversal 
cutoff dependence (in q z ) of the Feynman path integrals and their representation as a 
coherent state functional integral jlj . We interpret the integral over q z in such a way that 
its imaginary part vanishes. The real part is well-defined and simply gives a constant 
renormalization of the chemical potential. After integrating over k z and evaluating 
in the zero-frequency, long wavelength regime we find 



E(q Q)- A [™± k ' + ^' d > q - 

x ' 4 J (2n) d 2 ^KJ^\k ± \(^/KjK~ 3 \k ± \+Dkl-fi) 

A 2 fd d k ± 4(1 -l/% 2 [d d kj 



4 J (2n) d 2y/K^\k x \{y/K^/K^\k x \ + Dk\ - fi) J (2tt) 



(6.7) 



Although the correction to ft in (|6.5|) is well-behaved, the renormalization of D represented 
by the coefficient of g 2 diverges for d < 2, when Jh « 0. This infrared divergence suggests 
that renormalization group techniques be employed to study the system for d < 2. Our 
approach follows |48|], which is a variation of the dynamic renormalization group of 



First we integrate out a momentum shell in k± from Ae~ e to A, but integrate freely over k z . 
We then rescale our variables so that the ultraviolet cutoff is held fixed. After rescaling 
and 5n accordingly, we are left with the same theory but with different coupling constants. 
When this procedure is iterated, A and v are driven toward a fixed point which describes 
the universal long wavelength behavior in the dilute limit. 

6.1. Momentum Shell Integration 

We must first integrate out the transverse momentum in the range Ae~ l < q ± < A. 
This can be done straightforwardly by expanding the functional integral ( |6.1|) in v and A. 
The expansion can be represented diagrammatically as in Figure 12. Care must be taken 
to account for all possible contractions of the operators in the expansion. The symmetry 
factors can be found in the usual way for Wick expansions. It is important to note that 
diagrams renormalize the remaining low momentum modes, and are not simply expectation 
values. The diagrammatic rules may be summarized as follows 

• For each line, assign a momentum ki while conserving energy and momentum at the 
vertices. 
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For each polymer line include a factor of 
mined by the direction of the line. 



For each nematic line include a factor of 



-ik z + Dk\ — jJL 



. The sign of k z is deter- 



+ 



S ij - k l ± k J ± /k 2 ± 



K x k\ + K z kl K 2 k 2 ± +K 3 k 2 z 

• If a vertex joins four polymer lines, include a factor of 2v. 

• If a vertex joins two polymer lines and one nematic line, include a factor of ^ and the 
sum of the incoming and outgoing polymer momentum. 

• Divide by the symmetry factor. The symmetry factor is the product of the number of 
ways of permuting the vertices and the number of ways of permuting the lines while 
keeping the contractions the same. 

• Integrate over all momenta q z , but only integrating the transverse momenta from Ae~ £ 
to A. 

Upon carrying out this procedure, we arrive at the following relations for the inter- 
mediate values of the coupling constants: 



D 

A 

v 



D 1 + 



A 



2 r 



D 



K 1 (d-1) + K 2 



K\K 2 



A d A d ~ 2 (1 - e-( d -W) 



2d 



d-2 



= A 



v — 



X 2 v 



+ 



A 4 \ A d (l - e-( d -W) 



K 
/' - : /' 



AD AKiD ' 16KfD 
A 2 A A,DA d 



d-2 



i = 1,2,3 



v 

D + 8K!D 



d 



(1 



-de 



) 



(6.8o) 
(6.86) 
(6.8c) 
(6.8d) 
(6.8e) 



where A d = 2/[r(|)(47r) d / 2 ]. In evaluating the integrals we have ignored terms of two 
types. If a term diverges as i — * oo in a smaller dimension than the most divergent term, 
it is irrelevant by power counting, further, there are terms which appear to diverge in 
a higher dimension, but they are all higher order in the external momenta. When we 
rescale we must rescale these momenta. Doing so will render these terms irrelevant in 
the renormalized and rescaled theory. The intermediate values of D, p,, and the Frank 
constants reflect all the contributions to them at one loop order, while the expressions for 
A and v represent only the most relevant contribution to their values. 

These intermediate coupling constants have not yet been rescaled. We rescale 
lengths by V = Le~ £ and times by T' = Te Jo 1<yi ^ M . We now set q' ± = q ± e £ and 



= q z eJo 1<yi ^ M , where the function ^(i) is to be determined. The dimension of ip(r, t) is 
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just (L d )" 1 / 2 . Note that even when we rescale, ifj(r : t) will rescale trivially because it has 
no time dimensions. ijj(r,t) has no anomalous dimension in our renormalization scheme, 
to leading order in e = 2 — d. 

When doing the first momentum shell integration, the coupling constants were inde- 
pendent of length scale. However, they then acquire a momentum dependence because we 
have absorbed the large momentum effects into them. The correct renormalized theory 
is a coupled set of integral equations where the coupling constants are taken to be scale 
dependent. An alternative, but equivalent approach is to integrate over a small momen- 
tum range where the coupling constants are approximately fixed and then repeat the entire 
calculation iteratively. This leads to the usual differential renormalization group equations. 

6.2. Recursion Relations 

We now choose an infinitesimal momentum shell e~ s , and take the limit 5 — > 0. This 
leads to differential renormalization group equations which can be integrated to produce 
the couplings appropriate for a cutoff q ± < Ae~ e . We use units such that A = 1 in the 
following. The differential recursion relations are 

dD(£) / A 2 KAd- 1) + K 2 A d \ 

■ - D -2 + 7+ — u ' (6.9a) 



d£ V D K i K 2 2d 

^=M-l+7) (6.%) 

dvU) ( , v A A 2 , \ A 4 „ . n . 

d + 7-wr;A d +——A d ) - ^o—, A d (6.9c) 



dl V. 2D u K X 4D U J Kf32D 

dK x (l) 



dl 
dK 2 (£) 

di 
dK 3 (£) 

dl 



Kx (d - 2 + 7) (6.9d) 
K 2 {d-2 + 1 ) (6.9e) 
K s (d-"/) (6.9/) 



dm -^~(^^)A d D (6.9,) 



d£ n ' \D SKxD^ 

It is convenient to chose j(£) so that the renormalized, rescaled D remains fixed at 
its initial value. By examining ( |6.9j ) we can see that the there are three dimensionless 
coupling constants which come into this theory: Ai = A(ifi(d - 1) + K 2 ) 1/2 (K 1 K 2 D)- 1 / 2 , 
A 2 = A(Ki J D)" 1 / 2 and v = v/D. 

The quantity, Ai , controls the coupling between the tangent field of the polymers 
and the nematic matrix. The coupling, A2 , represents the interactions between density 
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fluctuations in the polymers and the the nematic. Note that none of the dimensionless 
couplings depend on the bend Frank constant, K%. For non-interacting polymers, q z ~ q\. 
Thus we expect that the term ^Ksq 2 \5n\ 2 ~ q A ± and will be suppressed at very long 
wavelengths. This is why the bend elastic constant does not couple to the theory at this 
order, though presumably it will at higher order. 

The result that the geometric mean of the splay and twist Frank constants comes into 
Ai is not unexpected. In || it is shown that the effective rigidity of a polymer in a nematic 
is the harmonic mean of its original rigidity and the induced rigidity from the nematic. 
We can understand the meaning of this harmonic mean by considering the following static 
example of a single polymer. If the polymer deviates away from the z direction, the nematic 
can relieve the stress in the same time slice by a twist, a splay or by a combination of the 
two. In d dimensions, there are d — 1 twist directions and only 1 splay direction. In 
the static limit, where q z = 0, the nematic bends back to its preferred direction over the 
transverse plane, whether the relief is through a twist or a bend. This allows us to estimate 
the nematic energy cost in disturbing the polymer, 

1 — IT 

5F oc a 2 K x + (d- 1)(——) 2 K 2 (6.10) 

d — 1 

where < a < 1 measures the fraction of the distortion relieved by a splay mode, and 
(1 — a)/ {d — 1) is the fraction of the distortion carried by each of the equivalent twist 
modes. Minimizing the energy with respect to a, we find that 

showing that the coupling to tangent fluctuations is given by A\. 

Finally, A2 is the coupling between the polymer and the splay degrees of freedom of 
the nematic. As we have seen, this coupling comes about through density fluctuations of 
the polymers. The effect of the nematic is to create an attraction between polymers, not 
unlike the Van der Waals attraction between neutral bodies. 

6.3. Fixed Points and Flows 

Our theory is described by the full space of all the coupling constants. We can describe 
the flow of the theory towards a stable theory by analyzing its fixed point structure in this 
space. The behavior of the Frank constants is trivial, amounting to a mere rescaling. We 
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assume, moreover, that the chemical potential has been adjusted to the critical point which 
describes the critical point of the theory 

We will analyze the flow in terms of the variables described above, namely v , A2 and 
Ai near d = 2. From ( |6.9| ) we find: 

4 



dv 



~d£ 

d\ 2 
dt 



A, 



64tt 



A, I 4 



V 

16tt 

A^ 
16tt 



Note that A2 is slaved to Ai in the sense that their ratio is independent of 

— 2 — 2 

u = v — (1/4)A2 , the recursion relation for u{£) depends only on Ai . 

2 



du 
d£ 



u e 



V 

8tt 



(6.12a) 

(6.126) 

(6.12c) 
If we let 

(6.13) 



— 2 — 2 

so that it suffices to consider flows in the space of u{£) and Ai {£). The subspace Ai = 

is the theory considered by Nelson and Seung for flux lines, while the subspace u = 

is the non-interacting theory considered by de Gennes [§. By examining the flow to the 

fixed point, we can decide which theory dominates in the long wavelength limit. 

Since the dimension of interest is d = 2, we discuss the flows for e = 0. In this case 

v 2, 



we can solve (16~T26D for X 1 (£) 



Ai (£) 



Ai 2 (0) 



8tt 



1 + 



Ai>) 
8tt 



(6.14) 



and then solve for u, 



u{£) 



u(0) 



1 + 



Ai 2 (0 ) 
8tt 




4n 
£\n£ 



(6.15) 



if Ai(0) 7^ 0. If Ai(0) = 0, u(£) is given asymptotically by 4tt/£, in agreement with the 

results in [J]. The flows are illustrated in Figure 13. 

- 2 - 2 

If wo > 0, the flows go the origin in the (Ai , •u)-plane. For nonzero Ai (0), they come 

- 2 

into the origin tangent to the Ai axis. This means that the logarithms that de Gennes 
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discussed for a single polymer dominate the logarithms associated with interpolymer in- 
teractions. If, on the other hand, uq < 0, the flows will still run to Ai = 0, but u will run 
off to — oo. Since u is the coefficient of p 2 in the hydrodynamic theory, this flow will push 
the system through a gas-liquid phase transition, since presumably there are higher or- 
der terms in both the hydrodynamic and boson language, corresponding to many-polymer 
interactions. 

Finally, we can calculate the logarithmic corrections to the wandering of a single 
polymer in the dilute limit, 

l\ (T\ rmia\ /^r 2 (^(r,L>/,*(0,0) ) 

r (L -r 0) = r\,*(n n\\ ■ 6,16 

J dr(V(r, L W (0,0)) 

In order to do this, we observe that 

( ^(r, z)V*(0, 0) ) e=Q = e M (iP(e- l r, e~ 7( ^'^*(0, 0) (6.17) 

We have chosen j(£) so that .D will remain fixed at every scale, i.e., 7(£) = 2 — Ai (£)/87T. 
This choice is arbitrary and, of course, does not affect the expressions for physical quanti- 
ties. 

We now choose such that z' = e J« z = ao, where ao is the persistence 

length of a single, directed polymer. As z — > oo, £* becomes large enough so that we 
have flowed very close to the stable fixed point in Figure 13. Since Ai, A2 and v are small 
near this fixed point, we can simply use the zeroth order term in perturbation theory to 
calculate the two-point function. We have 

( \t(z) - r(0)| 2 } e=Q = e 2r ( \r'(z') - r'(0)| 2 ) e=e . 

ar /dV(r0^(r',20f(O ) 0)) w . 



JdV(^(r',^*(0,0))^ 
2dDe 2r z' = 2dDzJo^ 2 -^ de ' 



(6.18) 



T 2 



Substituting Ai (£)/8n for 2 — 7(£), we can integrate (|6.14|) when d = 2 and find 



|rW-r(0)| 2 ) fco =4iJ Z 



8tt 



(6.19) 



Our choice of £* amounts to choosing £* w iln(z/ao)- Substituting this into ( |6.19[ ) and 
writing Ai in terms of the original couplings (with A = 1, as it was originally) 

( \t(z) - r(0)| 2 ) = ADz + f)l K K \ z^ (-) (6.20) 

47i(K 1 K 2 ) \a J 
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again resulting in a logarithmic correction to wandering. Comparison with de Gennes' 
result (|2.7| ) shows that the logarithmic correction to wandering is only half as large as 
would be predicted by the simple argument in Section 2. 

The correlation ( ?/?(r, z)if)*(0, 0) ) is the probability distribution for the wandering of 
a single polymer only in the dilute limit. This correlation function represents inserting 
a polymer at (0, 0) and removing a polymer at (r, z) , but there is no constraint that 
it be the same polymer. However, if the system is sufficiently dilute, the likelihood of 
polymers swapping their heads and tails is small. In the above derivation, we halted our 
renormalization group iteration when z ~ ao- However, we can also stop iterating when 
the polymer density, po, becomes of the order of the A 2 . Since A -1 is of order the monomer 
thickness, we can then apply the hydrodynamic theory of Section 5. 

For fixed z, we can always make the system dilute enough so that as we follow the 
renormalization group trajectory z(£) ~ ao before p(£) ~ A 2 , where z(£) = ze Jo 1{ ' % 
and p(£) = e 2i po. In this regime the wandering is given by ( |6.20|) . However, as the density 
increases, we come to a point where p{£) ~ A 2 before z(£) ~ ao- Now we must choose £* 
so that p(£*)A 2 = 1. In this case ( |6.20| ) will cross over to 

<| rW - r (0)|»> = 4D, + J^^ln^) (6.21) 

This result will hold when ^ > ^ 2 . In this very long polymer regime, the interpolymer 
interactions destroy the logarithmic correction to wandering. Each time the polymers 
wander into each other, the random walk is reset, and thus the logarithm does not build 
up along their length. 



7. Effects of Hairpins 

In the preceding sections we considered polymers without hairpin configurations. As 



in Section 2.3, we can account for hairpins by adding a term to the action (|6.2|) , namely 



S -+ S - I j dzd 2 r[ip 2 + (rf] (7.1) 



The coupling w oc exp(— e^/fc B T), where the hairpin energy is related to the coupling 
constants in ( p..2| ) by = O(^fgR) - see Section 2.3. Upon repeating the analysis of 
Section 3.2, we see that these terms create and destroy pairs of polymer lines and so add 
hairpins to the theory. 
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Upon carrying out perturbation theory in v and A, we find that unphysical diagrams 
appear in our theory. By allowing hairpins, we now also include loops of interacting 



loops of polymers by replicating ip(r, z) M times, and then taking the limit M — > 0. This 
could also have been done in our earlier analysis, although it is unnecessary, due to the 
retarded nature of the propagators. In the theory without hairpins but with free ends, 
one might think that the free ends generate an effective hairpin term as in Figure 14b. 
However, the effective hairpin strength is 0(h 2 ), and makes good physical sense. With free 
ends present, a long polymer may interact with two short polymers, simulating the effect of 
an intervening closed loop polymer. We must be certain that taking M — > preserves these 
graphs. Consider the action for M polymer fields coupled to a source (without explicit 
hairpins) , 



By changing the replica basis, we can choose to only allow the a = 1 direction to have 
a source term. Because of the preferential status of this direction, the "effective" closed 
loops will be 0(1), as opposed to O(M) in the case of real closed loops. Thus taking 
M — > will not alter the results found earlier for free ends. 

In the dense phase, hairpins lead to effects very similar to those discussed for free 
ends in Section 4. To see how hairpins affect the transition near Jx = 0, We first set 
tp = (tpi + iip 2 )/V2, and note that the action takes the form 



polymers as in Figure 14a. However, following de Gennes , we can eliminate closed 




> + 



F n [5n] 



(7.2) 



S 



F n [8n] 



iip a ±d z ip a2 + -4>ai[-DV 2 j_ - (fx - w)]ip, 



I 



(7.3) 



dzd 2 r I + -ip a2 [-DVl - (A + w)]ip a2 + -(iPl 1 + ip 2 M%+ip%) > 

+ i^-5n ■ (tp al V ± tp a2 - V^VxVai) 
v 2 / 
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where we have followed the usual summation convention. Recursion relations can be 
constructed as before for finite M, resulting in 



dfl(£) 
dl 

dw{£) 
dl 



u 



D 



(7.4a) 
(7.46) 



where u = v — A2 /4 was defined in the previous section. The recursion relations for the 
other variables remain the same as in ( |6.9a — /| ) . Upon taking M — > we see that only 
the recursion relation for p is different from fl6.9g ), and that (since 7 = 2) w is a strongly 
relevant perturbation. Although ft has a different recursion relation, this has no effect 
on the universal long wavelength properties. The M — > limit amounts to removing the 
second graph in Figure 7c, while keeping the first graph. As before, the surviving graph 
can be absorbed into a change in the chemical potential. 

Suppose our system is very dilute, i.e. ft<0 and w > 0, then %j) 2 will condense. We 
can follow ft + w via the recursion relations ( |7.4|) . The coupling ft + w grows rapidly under 
iteration, until it becomes large enough so that we can integrate out the massive modes 
corresponding to the {i/j a2 }. The resulting theory involves only ip at and 5n, 



( 1 



dzd 2 r 



+ y(VWa0 2 



2 2 ' T (V r n: -dn)- 



k + A(V^^ Q1 ■ dn)d z i(j al 



> + 



F n [5n] 



(7.5) 



with ft 1 = ft — w. We now have a theory in which directed propagators are replaced by 
merely anisotropic gradient couplings - the meandering in the parallel direction scales 
just as the meandering in the the perpendicular direction with different proportionality 
constants. The couplings shown in (|7.5|) only serve as a caricature of those found in 
the theory - they will be changed by numerical factors depending on at which point we 
integrate out ip 2 - As we go through the dilute-dense phase transition by letting fi 1 change 
sign, the theory describes a quasi-isotropic polymer melt with self avoidance in the limit 

m -> igg. 

The dilute limit of directed polymer melts without hairpins is described by an XY-like 
critical point (with a diffusive propagator), as indicated by Figure 8. The change in this 
phase transition induced by hairpins can be summarized by considering the mean field 
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diagram (see Figure 15) associated with the polynomial part of ( [7.3|) . We imagine holding 
v and w fixed, and varying the chemical potential ft in the plane of r\ = —p + w and 
r 2 = —p — w. When w = 0, the critical point, at r-y = r 2 = is just that considered 
in Figure 8. When w ^ 0, however, the trajectory as p varies passes through one of the 
Ising-like critical lines with a quasi-isotropic effective propagator (see ([7.5|) ), on the r\ — 
and r 2 — axes. In a similar fashion, the line describing the XY-like symmetry of the dense 
phase for t\ = r 2 < becomes unstable to Ising-like dense phases when w = 0. 

The above discussion ignores the coupling to the nematic field. This theory is similar 
to the compressible Ising model [pi] though the marginal operator present in that theory 
is not present in ours. This is guaranteed by the underlying nematic symmetry of our 
theory. Note that the original field theory is invariant under 5n — ► — Sn, ift *->■ ip* and 
z — ► —z. This prevents a term such as (V x -^n)-!/^, the usual coupling of the Ising model 
to an underlying elastic lattice. The couplings which do appear do not affect the above 
arguments, and we are left with an M component Ising model as M — > 0, reproducing 
de Gennes' theory of isotropic polymers. Note that the upper critical dimension of this 
Ising-like transition is d c = 4, as opposed to the result d c = 2 + 1 = 3 appropriate when 
u; = 0. 
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Appendix A. Consequences of Rotational Invariance 

The full rotationally invariant coupling between the nematic and the polymer is given 

by 

F = gjds 
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1 - 



dR(s) 
ds 



n 



+ F n [8n] 



(A.l) 



where s is the arc-length along the polymer. We now consider the following transformations 
of our fields (leaving all other fields the same), 



Tl x ► Tl x -\- V x 

n y -> n y + v y , (A. 2) 

R z (s) -> R z (s) - v x R x (s) - VyRy(s) 

where v = (v x ,v y ) is a constant vector in the plane perpendicular to z. It is easy to check 
that (|A.1|) is invariant under these transformations if H, the magnetic field, is 0. In terms 



of the small deviations defined in Section 2, this transformation, to linear order in v and 
the fields is 

8n'(r'(z'),z') = 8n(r(z),z) + v (A.3) 



and 

dr'(z') _ dr(z) 
dz' dz 

and thus to linear order in v 



(A.4) 



F' = F + g f dzv[^--8n{v{z),z)^ (A.5) 

However, we may absorb the shift in 5n into a redefinition of the field, so as to eliminate 
the effect of the transformation and find that F is invariant under ( |A.2| ) . 

Returning now to the coherent state field theory (|3.8| ) , we find that it must be invariant 
under 

^n^r', z') = 5n(r, z) + v 

(A.6) 

d' z =d z -v-V ± 

We now change the coefficient of ip*d z ift from 1 to some constant value a. The transfor- 
mations lead us to 



= S a [V>*,V,H + v- J dzd 2 r{-aWV^+ 1 -{r^^-^^) (A.7) 
= S a [il}*,il},8n} + (l-a) J dzd 2 rv -fV^. 

Thus S will only be independent of v if a = 1. 
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Appendix B. Derivation of Hydrodynamics 

In fact, the hydro dynamic theories described in Section 4 can be derived from the 
more microscopic boson theory ||52|| . To this end, we begin with (|4.3|) , taking h = to 
start, 



dzd\ 



D ( V±P ^ 2 + Dp(V ± 9) 2 + ipd z 9 + ipVJ -6n-fip+ -p 2 
4p ' 2 



(B.1) 



We now introduce a vector field P, via a Hubbard-Stratonovich transformation, in order 
to eliminate the Dp(V ± 9) 2 term in ( |B.1| ). We now have 



VPVpV9e~ s 'v 



with 



SL 



dzd r 



Ap 



+ ipd z + ipV±6 ■ Sn 



-pp+ -p 2 + DpP 2 + 2iDpP ■ VJ 
2 



(B.2) 



(B.3) 



If we integrate out P we return to the original action S p . However, we can now integrate 
over 9. Since it appears but linearly in S' p , its integration results in a delta-functional 



Z = / VPVpe- SH d[d z p + V ± ■ {pSn + 2DP) 



where 



S 



H 



dzd 2 r 



DpP z + 



2 , D(v ±P y 



Ap 



(B.4) 



(B.5) 



Thus we have traded interactions between 9 and 5n for a constraint relating p, P and Sn. 
We must now identify the physical meaning of P. In the case Sn = 0, the constraint in 
( p.4j ) becomes 

d z p + V ± ■ (2DpP) = (B.6) 

which suggests that 2DpP is the "tangent" field introduced in Section 5. However with 
the nematic field included, we would have 



d z p + V x • (2DpP + pSn) = 0. 



(B.7) 



Returning to the path integral for a single polymer ( |2.1|) , we identify the Euclidean La- 
grangian as 
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Leading to the Euclidean momentum 



dL i ( dr ^ , 



dr 2D V dz 
so that 

— =2Dip + Sn. (B.10) 

az 

Apparently, P is the field associated with ip. Letting v be the field equivalent of the 
velocity, we have 

v = 2DP + 8n, (B.ll) 

so our constraint becomes 

d zP + V ± -( P v) = (B.12) 

Note that pv = t is the "momentum" field of our quantum bosons, which we have called 
the "tangent" field in the previous sections. We now replace P by (v — 5n)/2D in ( |B.5| ). 
Upon expanding p around some mean field value po, we are lead to the hydrodynamics of 
Section 5. 

If we were to add sources to the theory, ( P-l|) is changed by 



5S P = J dzd\2h^fp~ cos (B.13) 

We may expand this term in powers of 6 if the fluctuations in 6 are small, which they will 
be if the polymers are sufficiently long. To lowest order ( [B.13 ) is 



SS P = j ' dzd^r h^/p-Q 2 (B.14) 
Now upon integrating out 9, we no longer have a 6 functional as in ( |B.4| ), but instead 

VvVpe~ Sn (B.15) 



with 

D „(v-^ + ^!- P „+v 

S„ = I dzd\ { P ) (B.16) 

resulting in the hydrodynamic theory of finite length polymers ( |5.5| ). There we had the 
term 

\^T [dzP + V ' t)2 - (R17) 
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Recalling G = k B T£/2p and I = p 1//2 /h, we have ^(G/k B T) = (4 v /po/i) _1 , as required. 

The integration over 9 in ( |B.3| ) needs a more careful analysis. Indeed, it will constrain 
the number of lines passing through a surface to be integer. When changing variables 
from if; and ip* to p and 9, theta is only determined modulo 2n. Because of this, we can 
consider field configurations in 9 which increase by 2n as we go around a particular line. 
This vortex line is analogous to a point vortex in the two-dimensional XY model. We 
can now consider a closed vortex loop (the three dimensional analog of a vortex-antivortex 
pair) where 9 jumps from to 27m as we traverse a surface with the vortex loop as the 
boundary (a so-called "branch" disc). Rewriting the 9 dependent part of (p3.3|) , we have, 
for a volume E bounded by the surface 9E, 

z / dV (pd z 9 + t- V ± 9) =i [ dS-fd-i I dV 9 (d z p + X7 ± ■ t) (B.18) 

where T = (t x ,t y ,p) is the three-dimensional "particle" current. The volume integral on 
the right hand side of (|B.18| ) will lead to the constraint (|B.12| ) in the volume V, when 9 is 



functionally integrated over. The surface integral contains the new constraint. 

Consider a vortex line which lies on the surface <9£ (see Figure 16). If we now split 
<9£ into the parts above and below the surface of discontinuity, dY* u and <9£d, respectively, 
the surface integral is 



2iiin / dS-T + i dS-T9 s (B.19) 

where 9 S is the smoothly varying part of 9. In the same way that we generated the equation 
of continuity, the second integral will constrain the total flux through the closed surface 
<9£ to be (no free ends.). However, upon summing over all possible values of n, the first 
integral in (|B.19|) will constrain the "flux" of lines passing through the surface <9£ u to be 
an integer. If polymers are dense and entangled, this set of nonlocal constraints should be 
unimportant in the hydrodynamic limit, where the fluctuations in 9 are small. 
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System 


Mechanism For 
Alignment 


Free Ends 
Important? 


Hairpins 
Important? 


Polymer 
Nematics 


Spontaneously 
Broken Symmetry 


Yes 


Yes 


Polymers In A 
Nematic Solvent 


Induced By 
Solvent 


Yes 


Yes 


Ferro- and 

Electrorheological 

Fluids 


External 
Fields 


Yes 


No 


Flux Lines In 
High Temperature 
Superconductors 


External 
Fields 


No 


No 



Table I. Characteristics of the different physical systems 
considered in this this paper. 
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Figure Captions 



Fig. 1. Hydrodynamic averaging volume surrounding a small region of a polymer nematic 
containing many polymer strands. The average over the polymer tangents in this 
volume defines a coarse-grained director field n(r,z), which then tends to align 
the individual polymers which pass through the region. 

Fig. 2. Conventional short chain nematogens (ellipsoids) connected by hydrocarbon spac- 
ers to make a polymer nematic. 

Fig. 3. Isolated polymer in a short chain nematic solvent. The Goldstone modes of the 
nematic produce anomalous wandering of the polymer transverse to the z-axis. 

Fig. 4. Polymer interacting with director fluctuations in a two-dimensional nematic 
medium, and described by its complex coordinate z as a function of arclength 
s. 

Fig. 5. A directed polymer propagating either up or down the z-axis, changing direction 
whenever a hairpin is present. Each hairpin contributes a factor of w/2, and each 
line contributes a free propagator, Go to the full propagator G. 

Fig. 6. Contribution to order h 8 to the polymer generating function. Solid lines are the 
polymer propagators Go(q ± , q z ). Dashed lines represent the interaction potential 
V(r), while the dotted lines represent interactions induced by the background 
nematic field. 

Fig. 7. Ingredients of the graphical perturbation theory described in the text. Figures 
(a) and (b) represent contributions to the renormalized interpolymer interaction 
potential which vanish due to the retarded nature of the polymer propagator. 
The terms shown in (c) are constants which can be absorbed into a shift in the 
chemical potential. 

Fig. 8. Phase diagram as a function of the boson order parameter and the chemical 
potential p,. The dashed line represents a contour of constant polymer length. 
Typical configurations are shown at the points A (dilute) and B (dense and 
entangled). The solid curve represents the boson order parameter for h = 0, and 
terminates at a critical point which describes the theory in the dilute limit. 

Fig. 9. Contours of constant scattering intensity for long chain ferro- and electrorheo- 
logical fluids. The contours are linear near the origin and surround a maximum 
on the q^-axis located approximately at the position of the first Bragg peak of 
the nearby hexagonal crystalline phase. The linear contours near the origin are 
rounded by finite chain length effects. 
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Fig. 10. Contours of constant scattering intensity near the origin for polymer nematics. 

The characteristic square root contours are rounded off as indicated by the dashed 
lines unless the polymers are very long. 

Fig. 11. Lowest energy solid phase contribution to the correlation function G(r, r'; z, z'), 
which inserts a flux head and tail into a crystalline vortex array. Dashed lines 
represent a row of vortices slightly behind the plane of the page. In (a), a va- 
cancy is created at "time" z, which then propagates and is destroyed at time z' . 
Interstitial propagation from z' to z is shown in (b) . The energy of the "string" 
defect connecting the head to the tail increases linearly with the separation in 
both cases and leads to the exponential decay of G(r, r'; z, z'). 

Fig. 12. Lowest order contributions to the perturbation expansion used in calculating 
corrections to the zeroth-order parameters in the theory. Graphs which are iden- 
tically due to the retarded nature of the polymer propagator are not included. 
In this case, the internal lines are only integrated over a small spatial momentum 
shell e~ l < q ± < 1, but are integrated over all "time" -like momenta q z . Figures 
(a) and (b) shows the contributions to the self-energy of the polymer propagator. 
The contributions in (b) are just constants which we must absorb into a redefi- 
nition of the chemical potential, £L. In (c) we show the graphs which renormalize 
the four-point coupling v. Figure (d) shows the corrections to the vertex between 
the nematic director 5n and the polymers. Figure (e) shows two graphs which 
contribute to the four-point coupling but do not vanish identically. However, 
their contribution is irrelevant by power counting. 

Fig. 13. Fixed point flow in 3 dimensions (d=2). We show u, the effective hard core 

- 2 - 2 

repulsion and Ai . The flows come into the origin tangent to the Ai axis. If 
uq < then the flow in uq runs off to large negative values, and the system will 
go through a gas-liquid phase transition. 

Fig. 14. In these graphs, we represent a hairpin insertion by a solid square and the insertion 
of a free end by a solid circle. Figure (a) shows diagrams which involve closed 
loops of polymers interacting with physical polymers. These graphs do not factor 
into a physical part and an unphysical part. Figure (b) shows graphs appearing 
in a theory with free ends which appear to involve closed loops of polymers, but 
instead are interacting with two or more short polymers. These graphs are indeed 
physical. 

Fig. 15. Trajectory with varying ft through the r\-r<i plane, where r± = —ft + w and 
T2 = —ft — w. For w 7^ 0, the transition from the dilute limit corresponds to an 
Ising-like phase transition. The line r\ = r2 < corresponds to an XY-like phase. 
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Fig. 16. Surface on which a vortex loop lies. The plane, denoted by hashed lines, is the 
"branch disc" on which the phase jumps discretely by a multiple of 2n. 
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